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Abstract. Let r be a cofinite Fuchsian group acting on hyperbolic two-space H. 
LetAf = r\H be the corresponding quotient space. For7,acIosed geodesic of A/, 
let Z(7) denote its length. The prime geodesic counting function 7rj\/(it) is defined 
as the number of F-inconjugate, primitive, closed geodesies 7 such that e'(T' < u. 
The prime geodesic theorem implies: 

/^3/4 " 

0<Am,j<1/4 

where = Aj\/_o < -^A/,i < ■ ■ ■ are the eigenvalues of the hyperbolic Laplacian 
acting on the space of smooth functions on M and = 5 + \ ~ ^M,j ■ Let 
Cj\/ be smallest implied constant so that 



„3/4 

< Cm for all u> 1. 

log u 



0<Aju,j<l/4 

We call the (absolute) constant Cj\/ the Huber constant. 

The objective of this paper is to give an effectively computable upper bound of 
Ci\/ for an arbitrary cofinite Fuchsian group. As a corollary we estimate the Huber 
constant for PSL(2,Z), we obtain Ca/ < 16,607,349,020,658 » cxp{30.44086643). 



Introduction 

Let r be a cofinite Fuchsian groupQ, and let M = F \ H be the corresponding 
hyperboHc orbifold. Let C{M) denote the set of closed geodesies of M, and let 
V{M) denote the set of prime (or primitive) closed geodesies (see IBus92i page 
245]). For each 7 € C{M) there exists a unique prime geodesic 70 and a unique 
exponent m > 1 so that 7 = 7™. Let /(7) denote the length of 7. Associated to 7 is 
a unique hyperbolic conjugacy class {^'7}r with norm 



The prime geodesic counting function ttm{u) is defined to be the number of F- 
inconjugate, primitive, hyperbolic elements 7 £ F such that e'^ < u. 

Selberg IISel56l and Huber IIFiub59[ IFiub61a[ IHub61bl independently proved the 
prime geodesic theorem, which is the asymptotic formula 

7rAf(u) ^ r-^, {u-^co). 
log It 
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Later, Huber proved a stronger version of the prime geodesic theorem, with error 
term^ 

0<Am,3<1/4 ^ ^ ■ 

where = Aj\f^o < ^m,i < ■ ■ ■ are the eigenvalues of the hyperbolic Laplacian 
acting on the space of smooth functions on M and 



For a proof of this theorem, via the Selberg zeta function, see |Hej76 Hej83| and 



IIRan78l . Randol IIRan84l and Sarnak llSarSOl also gave of proof using the Selberg 
trace formula. We call the implied constant inherent in Om{^^) the Huber con- 
stant, denoted by Cm . The Huber constant is an absolute constant, not an asymptotic 
constanfl The Huber constant is the minimal constant that satisfies 



0<AAf,j<.l/4 



<C'm{- 1 forallu>l. 

lORM 



The main goal of this paper is to estimate Cm for an arbitrary cofinite Fuchsian 
group. We determine an effective algorithm by which one can obtain explicit 
boimds for Cm using elementary geometric and spectral theoretic data associated 
to M. In particular, we apply the algorithm in the case of the full modular group 
PSL(2, Z) and obtain a precise, numerical bound for the Huber constant. 

In | |JK02| , the authors studied the Huber constant in the following situation: Let 
Fq be a cofinite Fuchsian group, and let F be a finite-index subgroup of Fq. They 
showed that 

Cm < [Tq : F] • C]\fg. 

Hence if one could estimate this constant for Fq = PSL(2, Z), one would have an 
estimate for any congruence subgroup. For more applications of these ideas, see 
| JK09l|jK06l|jKnTl . 



In the articles |JK09| , |JK06| and |JK01| , the authors studied analytic aspects of 
Arakelov theory, ultimately developing bounds for special values of Selberg's zeta 
functions. Green's functions and Faltings's delta function. The bounds for these 
analytic functions involved many explicitly computable analytic and geometric 
quantity, together with the Huber constant. As such, specific, effective bounds 
for the Huber constant then could be used to make the bounds derived in |JK09| , 
| |JK06] and [JKOIJ effective. With this said, the results in the present paper complete 
the analysis in the aforementioned articles, thus providing effective, numerically 
computable, bounds for various analytic quantities which appear in the Arakelov 
theory of algebraic curves. 

Main result. Our main result is Theorem 14. 151 (see ^A.8i , where we determine an 
upper bound for the Huber constant Cm in terms of various well studied invari- 
ants of F and M such as: the number of small eigenvalues A, the number of ex- 
ceptional poles of the scattering matrix J\f, the smallest positive eigenvalue Ai , the 
length of the smallest closed geodesic of M, the decomposition of the fundamental 



^Ruber's error term was slightly different: O (log it) ~^/^) 
'^If Cm was an asymptotic constc 
estimates for the error term Ilwa02l 



'^If Cm was an asymptotic constant, it would equal zero for the Modular group since there are better 
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Table 1 . Fundamental constants 



Constant 


Value 


i 


cofinite Fuchsian group with fundamental domain J~ 


T 


IlUIIlDcr OI lIlcClLllVcilcIlt UdrdUUllC CLlSpb 


1 


L , rill 1 ' ^ "i /-I / "V v"v> v~\ /"\ c ' 1 -1-1 / "\ v~\ /^-f 1 -v-i -1-/^ 1 c c 4-/^ r? X* i n /-I /-^/^t"v^ v~\ --i 4- X* , , 

— jy u lj^_Q 2 ' ueLoiiipobitioii or iiuo cu&p seciorb aiici Loiiipaci set j^y 


Or, iTiB, 






r^Ul t^Ut^UlllpclCL 1 Ci is Lilt: 1 1 ly Ucl UUlit^ LllcllllcLt;! (JL J 


A 


number of eigenvalues A of A, such that A < |, including Aq = 




1 2^"\/4 •^i ^i!i4 ]-^gj-g is the smallest non-zero eigenvalue 

i oo Ai > i, 


N 


(Ti > • • • > cTjv' are all the poles of 0(s) in the interval (15, 1] 


Cl 


minimum of the lower left hand matrix entry of non-parabolic elements 


c 


chosen to satisfy 1 < c < iV(r) for all hyperbolic 7 G F 


M 


log c, a positive lower bound for the length of all closed geodesies 



domain of F into cusp sectors, the area of and a few other simple, easily com- 
putable, invariant^ We could give the main result by stating our upper bound 
for Cm ■ However the result is so complicated that it would take us two pages just 
to list out the equation. So, we present our results as an algorithm, rather than a 
formula. To help the reader appreciate the complexity (for the case of a general 
cofinite Fuchsian group), we state, as a corollary, the simple case of a cocompact, 
torsion-free Fuchsian group. 

In the process of proving Theorem 14. 15[ we give an upper bound for the implied 
(absolute) constant of the spectral counting /i(r). See Theorem |3.3l In other words, 
for an arbitrary cofinite Fuchsian group, we find an explicit constant C, depending 
on r. so that 



/^(r) = |{r„|0<r„<r}| + 



47r 



-(2+^0 



dt<C 



(r > 0). 



Next, consider the following table of constants, which are valid only for F cocom- 
pact, torsion-free. 

Theorem. Let T be a cocompact, torsion-free Fuchsian group. Then the constant Cu is 
an upper bound for the Huber constant Cm, and the constant C is an upper bound for the 
implied constant of the spectral counting function of IS. . 

In order to estimate Cm for general cofinite Fuchsian group we need the following 
modified table. Some constants from the cocompact table are redefined. 

Theorem. Let T be a cofinite Fuchsian group. Then the constant Cu is an upper bound 
for the Huber constant Cm, cmd the constant C is an upper bound for the implied constant 
of the spectral counting function of A . 

To show the usefulness of Theorem 14.151 and Theorem 13. 3[ we estimate their re- 
spective constants (C and Cm) for the modular group PSL(2,Z). The results are 
Theorem l3.4l and Theorem l4.161 



These invariants are all known or easily estimated for the Modular Group. 
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Table 2. Algorithm for cocompact, torsion-free case 



Constant 


Value 


Ci 


2e - 2, e = 2.718... 


ID 


^^'^ (/lie 1 Q m oi"OT' A ^ 
ijpi \il lb Lilt: ClicilllcLfc;! Ul ^ ^ 


c 


3(m + 745s) 


A 


number of small eigenvalues of A including Aq 


Cio 


8480^/^ 


Cl2 


{A - 1)(1 + 3Ci + + Ci)) + 2Ci + 2 


Cl3 




Cl6 


C12 + Cl3 + ^ 1-7^1 ClO 


Cn 


4^ + 4Ci6 


^ io 


4^4, + 5Ci fi 


C 


1 < c < Ar(P) for all P e T, hyperbolic 


M 


logc 


Cl9 


•-^18+ 1-1/c 


^20 




C21 


|c-2|i5|2 + |2-V^ll5|7 


C22 


l-l/log{2) 


Cu 


C21»4 + C20j;;^ + C20 + (^20^^22 + |C'2oC'21 



Theorem. Let T — PSL(2, Z). Let ^(r) be the spectral counting function l|3.7b , and let 

C = 1,682,997. 

Then for r > 0, 

/i(r) < C (^r^ + . 
Theorem. Let T = PSL(2,Z). Then Cm < 16,607,349,020,658. 

Using the theory laid out in ISar82l , and the excellent open-source computational 
program PARl/GP BTheOSI , we can explicitly calculate (and list out) the length 
spectrum for PSL(2, Z). Our computations suggest that 

Cm < 2. 

So it seems that our result is not sharp. 

Why is our result for PSL(2, Z) so much larger than the experimental result? There 
are many reasons. Our bound for the spectral counting function is probably one 
million times too big (when compared to Weyl's law, an asymptotic result). We are 
mostly interested in how the Huber constant grows as a function of the various 
fundamental constants. Estimating the spectral counting function required abso- 
lute estimates of various infinite sums over primitive hyperbolic classes; and we 
used trivial bounds for the hyperbolic counting function. The non-trivial absolute 
bound being the sole objective of this paper. However, using PSL(2,Z) specific 
bounds on the spectral counting function, the result could be markedly improved. 

Connection to existing literature. As stated above, the boimd for the Huber con- 
stant obtained in the present paper allows one to complete the analysis in |JK09| , 
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Table 3. Algorithm for the general cofinite Fuchsian case 



Constant 


Value 


C^ 
^1 


2e- 


2 


r> 

JD 




Lemma[l.2lfor more details 


L-3 


8r 


Ci 


2m R sin Or 


{2e„ ^ 2(TT-eR,) ) 


C5 


8t + C4 + ^ + 745B 




C5+ (2|logci| + 


('^Ar-i)= j 


V5 


Cv 


,l(2|logc,| + ^) 


c 


3^6 4 


-Cj 


A 


number of small eigenvalues of A including Aq 




8480^ 




C12 


(^-l)(l + 3Ci + ^ 


-(l + Ci)) + 2Ci+2 


Ci3 




C'lo 


C 
Ul4 


^10~37r 1" 


§ + 2Tlog2 




56Cio 
3 






Cie 


C12 + Ci3 + Ci4 4- 


Ci5 + ^ 




Cn 


4A + 4Ci6 


^18 


iA + bCie 


C 


1 < c < N{P) for all P £ T, hyperbolic 




lege 


C'lg 


C18 + 


C20 




C21 


l^-2|ish + 




C22 


i 

l-l/log{2) 


Cu 


C21A + C2Qf^^ + (^20 


+ C'2oC'22 + JC20C21 



|JK06| [ and |JK01| , resulting in explicit and computable bounds for certain ana- 
lytic quantities in Arakelov theory. We refer the interested reader to these articles 
for the precise statements of the main results, from which one can immediately 
employ the bounds obtained here in order to further strengthen the analysis pre- 
sented in IJK09I , I IJK06I I and | |JK01| . We will discuss further applicability of the 
present paper to other computational problems. 

In | |ECJMB| , the authors used ideas from Arakelov theory and developed a strat- 
egy to compute Galois representations modulo £ associated to a fixed modular 
form of arbitrary weight. Ultimately, the goal of |ECJMB| is to devise an algo- 
rithm, which has complexity that is polynomial in i. A summary of the ideas from 
| |ECJMB| is given in IIEdi06l , where the purpose is to focus attention to the weight 
twelve modular form associated to PSL(2, Z). The ideas from ]ECJMB| have been 
used to achieve advances in many other computation problems; see, for example, 
llBoOTi, ICL091 , IICou09l , and iLaO^ for specific results, as well as the survey article 
llChlOSi for more general comments. 

As stated in ]ECJMB| , a key component of the algorithm was to determine bounds 
for Arakelov Green's functions. F. Merkl provided sufficient bounds to complete 
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the algorithm; however, the bounds from | |JK06| are sharper and, as a result, im- 
prove the efficiency of the algorithm from |ECJMB| . Theorem 4.8 of |JK06| sum- 
marizes a bound for the Arakelov Green's function in terms of hyperbolic data, 
including a special value of the Selberg zeta function. This special value was stud- 
ied in |JK09| and now can be bounded explicitly using the bounds for the Huber 
constant developed in the present paper. Thus, the bounds obtained for the Hu- 
ber constant proved here ultimately can be used in the algorithms constructed in 
I ECJMB| and subsequent articles. 

Going forward, it may be possible that the algorithms of | |ECJMB| , or similarly de- 
veloped ideas, could yield further applications, in which case effective bounds for 
the analytic aspects of Arakelov theory would play a role in the implementation 
of the resulting algorithms. Indeed, there has been recent progress in expand- 
ing Arakelov theory to allow for singular, or non-compact metrics (see IIBKK05L 
IIBKK07L IIFG09L and |Ha09l ). As a result, one can anticipate further compu- 
tational algorithms which utilize these recent extensions of Arakelov theory, in 
which case the bounds from |JK09| , |JK06 1 and | JKOl | and the present paper may 
play a role. 



Outline. In brief, the present paper amounts to a careful analysis of the proof 
of the prime geodesic theorem, making sure that each step in the proof yields 
explicitly computable constants. However, straightforward the methodology may 
sound, there are many details which involve considerable analytic difficulties. 

In ijl]we prove some useful lemmas and develop necessary notation. Our argu- 
ment consists of many estimates, each with an implicit constant, and our notation 
helps to keep track of these constants. At the heart of our argument is a variant of 
Karamata's Tauberian theorem (Lemma ll.ll in ^1.2|l that gives an absolute, instead 
of an asymptotic, result. 

In ||2]we study the heat kernel ^ * of the Laplace-Bel trami operator for the hy- 
perbolic orbifold A/ = F \ H. Our main tool is the Selberg trace formula. We 
specifically show that e~^^ = 0{l/t) for t e (0,5], and find the corresponding 
explicit 0-constant. Note that we focus on the interval (0,5] because small eigen- 
values (A < 1 /4) are treated differently in many applications of the Selberg trace 
formula. We will estimate the spectral coimting function ii{x) for A > 1/5 giving 
us a small margin of error from the 1/4 cutoff. 

Next, in f|3l using results from f|2]and our variant of Karamata's Tauberian theorem 
(Lemma ll.l|l , we prove Theorem 13.31 which provides an effective boimd for the 
spectral counting function. In i i3.2l we specialize Theorem 13.31 to the case of F = 
PSL(2,Z). 

Finally, in ||4l using Theorem l3.3l we obtain an upper boimd for the Huber constant 
(Theorem l4.15|l and in ^14. 91 we specialize to the case of F = PSL(2, Z). 



Concluding remarks. We conclude the introduction by addressing a number of 
points which naturally arise when studying the computations presented in this 
article. 

First, it is possible to further contract the algorithm for bounding the Huber con- 
stant and obtain a bound which can be related to the geometry of AI . For example, 
the constant ci is the radius of the largest isometric circle in a certain fundamental 
domain for M ; see page 38 of IIwa97ll . Also, the number of small eigenvalues can 
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be bounded by the topology of AI; see IZo82l . However, to employ general results 
would only further worsen the bound for Cm for special important cases. For ex- 
ample, for the congruence subgroups, the Selberg 1/4 conjecture asserts there is 
only one small eigenvalue, whereas the boimd from [iZo82J implies that the num- 
ber of small eigenvalues grows at most like the genus of the surface, which may 
be true in general cases but most likely is not true in special cases of particular 
importance. 

Second, the heat kernel style of test fxmctions which we analyze could be replaced 
by studying a different approach to the prime geodesic theorem, namely by fi- 
nite difference methods applied to the test fimctions denoted by tpk', see |Hej83 |. 



In our investigation, we saw no significant difference when looking ahead to the 
problem of bounding Cm, so we simply selected one approach. We will leave the 
problem of studying other ways of proving the prime geodesic theorem to see if 
the methodologies lead to improved bounds for the Huber constant. 

Finally, it is somewhat striking to see that in the case of the full modular group the 
bound we prove is approximately exp(30.45) whereas the numerical investigations 
lead one to believe that the true bound is log 2, which differ by an order of 10^^. In 
Fall 2003, the second named author (J.J.) was giving a talk at the Courant Institute 
at NYU on preliminary results associated to | JK01[ . (As described in |ECJMB| , it 



was then that Peter Sarnak communicated to the second named author (J.J.) the 
problem of bounding Arakelov Green's functions.) During the lecture, when the 
results of [JK09| were described, Sarnak stated that he felt it would be difficult to 
obtain a precise, numerical boimd for the Huber constant for any group, including 
the full modular group. With this stated, the present paper successfully addresses 
the problem of determining precise, effective bounds for Cm, and it remains to be 
seen how, or rather if, methods can be developed to refine the bound for the full 
modular group. 

1. Preliminary lemmas 

In this section we group together certain lemmas that will be used throughout this 
paper. Our key lemma is a variant of Karamata's Tauberian Theorem. 

1.1. Big 'O' notation. Let i E N, and let /, g be fimctions with a common domain 
D. The statement 

means that there exists a positive constant Ci so that 

\f{x)\<CM^)\ VxeD. 

If the domain in question is not clear, we will explicitly state it. We call the constant 
d absolut^ 

1.2. Karamata's Tauberian theorem. In this section we prove a variation Kara- 
mata's Tauberian Theorem llWid41[ Page 189-192]. Our version keeps track of the 
implied constants in the original theorem. 

We will apply the lemma below with the parameter d = 5. Small eigenvalues 
(A < 1 /4) are treated differently in many applications of the Selberg trace formula. 



^As opposed to an asymptotic constant, say Ai that would satisfy \f(x)\ < Ai\g{x)\ for x suffi- 
ciently large. 
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We will estimate the spectral counting function //(x) for A > 1/5 giving us a small 
margin of error from the 1/4 cutoff. 

Lemma 1.1. Let a{t) be a non-negative, non-decreasing function, continuous on all but 
a possibly countable subset of{0, oo). Suppose that a{0) = and that the integral 

POO 

f{s)= / e-''da{t) 
Jo 

converges for all s > 0. Suppose further that for C > 0, d> 

(1.1) /(s) < -, forallO<s<d. 

s 



Then 



a{t) < 3Ct for all t> i. 



Proof For x e [0, 1] let 

0<x<e-i 

1 e^^ < X <1. 



Note that 

g{x) < 3 (0 < a; < 1). 
Since a{t) is non-decreasing, we have, for s > 0, 

/•oo poo 

/ e-''g{e-'')da{t) < e-''3da{t). 
Jo Jo 

Next, observing that 

( 0<w<l 
) = | Kz., 

an elementary calculation using the fact that q;(0) = shows that 

e-''g{e-'') da{t) = / e-'^gie-"') da{t) = / e'^'e'' da{t) = a(l/s)-a(0) ^ a{l/s). 
Jo Jo 

Hence we have 

/•oo /"OO o*^ 

a{l/s)= e-'"g{e-'')da{t)< e~'' 3 da{t) = 3f{s) < — (0 < s < d). 
Jo Jo s 

Setting t = 1/s implies that 



a(t) <3Ct {t>-, 
' d 



□ 



Later on, we will let a{t) be the spectral counting function for the Laplacian A 
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1.3. Trivial counting lemma. Let F be a cofinite Fuchsian group. Let ■jt{x) denote 
the number of primitive conjugacy classes of hyperbolic elements (of F) of norm 
not exceeding x (the details are given below). In this section (following ||Fis871 
pages 47-49]) we give an explicit constant B so that for all x* > 1, Tr{x) < Bx. This 
result is called the trivial counting lemma^ 

Let p{z, w) denote the hyperbolic distance on H. For T hyperbolic, define the norm 
ofT, ^■(T),by 



N(T) = cxp inf p(z,Tz) 

Note that the norm is constant within a fixed conjugacy class {T}r- For .t > 1 
define 7r(a-) by 

7rix)^\{{T}reT\NiT)<x}\. 

Let T = TyHTi r\J^2 ^ ■ ■ • n J"^ be a fundamental domain for F, where each Tj' is 
a cusp sector and J^y is compact. Let Ci , • • ■ , Ct be a complete list of representatives 
of the F equivalence classes of cusps, and let Ai, . . . ,Ar be elements of PSL(2, M) 
satisfying Aj(j = oo. In addition, assume that AjFj = [0, 1) x [y, oo). 

If T = 1, then it follows that AiF has a positive Euclidean distance e to the real 
axis (e < Y). If r > 2, set 

a b \ . ._! 
c d j=^2A^ 

Note that c 7^ since A2A^^ does not fix 00; in this case set e = niin{Y, c~^Y~^}. 

Let Kj = {z E AjT \ e < Im(z) < F + 1}. It follows that Kj has compact closure, 
and hence has finite area \Kj\. Finally let dj be the hyperbolic diameter of Kj . 

Lemma 1.2. IFis87[ pages 47-49] Let The a cofinite, noncompact Fuchsian group. Then 
for all X > 1 

min(2,r) ^ 

- I] (cosli(loga; + 2rf,) - 1) 

Note that tt{x) = for x < 1. 
Lemma 1.3. Let 

min(2,r) 

B ^ At:Y{Y + 1) e^"^'. 

i=i 

Then 

tt{x) < Bx {x> 0). 

Proof. Since tt{x) = for x < 1, we assume that x > 1. By Lemma ll. 21 we must 
estimate 

min(2,T) 

T^(cosh(loga; + 2dj)-l). 

Note that for each j, [0,1) x [Y,Y + I] a Kj. Hence 

\KA > ^ 



^' - Y{Y + 1)' 



^Note that the hyperbolic (or geodesic) prime number theorem implies that ■k(x) = o{x) as a; — >• 00. 
Hence this lemma is called trivial, the non-trivial being the prime geodesic theorem. As our goal is an 
absolute constant, the trivial lemma will be used in a bootstrapping argument. 
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The lemma now follows by observing that for y, d > 0, 

cosh(y + 2d)-l< cosh(?/ + 2d) < ^(6^+^^ + 6^+2^^) = 6^+^'*, 
and by setting y = log x. □ 

1.4. Selberg trace formula. The main tool for this paper is the Selberg trace for- 
mula (see |Hej76|[Hi]83| , IIwa02l . IIVen82l ). 

Theorem 1.4 (Selberg trace formula). 

^ 1 f (p^ 1 1*^1 f 

h(rn) + — / h(r) — —(- + ir)dr—- — / h{r) r tanh(nr) dr 
f-' 47r 7r 2 47r 



n=0 

2re 



r 1 ^(7)^ -^(7) ' fTt 2TOfl sm 6'ij, 7r 1 + e '^'^'^ 4 2 



- T5(0) log 2 - — / h{r)%jj{l + ir) dr. 
Here A„ = i + are the eigenvalues of the hyperbolic Laplacian A and 



+ 1 



\-K if < A„ < \ 



Y An — -J if A„ > -J 

(see fJUfor some background material on A); h{z) is an even holomorphic function 
in the strip | Imz| < 5 + e satisfying 

\h{z)\^0{{\z\ + l)-^-y, 

{7}r are the conjugacy classes for the hyperbolic elements of F with primitive 
generator 70; (/)(s) is the determinant of the scattering matrix (3(s); is the hy- 
perbolic area of the fundamental domain F of F; g{r) is the Fourier transform of 
h{z); {R]r are conjugacy classes for the elliptic elements of F with 771^ = |C(i?)| 
(the order of the centralizer of i? in F) and tr(i?) — 2cos(0fl), < 0/? < tt; t is 
the number of cusps of F; and ^(s) is the logarithmic derivative of the Gamma 
function. See see )Hej76 Hej83| , IIIwa02l , IVen82l , for more details. 

2. Heat kernel for the hyperbolic Laplacian 

Let F be a cofinite Fuchsian group, and let H(F) denote the Hilbert space of mea- 
surable functions / : H 1-^ C satisfying 

• f[^z) ^ f{z) for all z e H, 7 e F 

• Xf 1/(^)1 dv{z) < 00. 

Let A : 'H(F) i-^. 'H(F) be the self-adjoint extension of the (two-dimensional hyper- 
bolic) Laplacian. In this section we will estimate (for < t < 5) 

For t > 0, let 



47r /r 6 2 

n=0 ^ 



h{r) = e 



-tr-' 
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9{x) 



Plugging in the pair /i, g into the Selberg trace formula gives an explicit formula 
for the (regularized) trace of the heat kernel * : 



— j e r tanh(7rr) dr 



logNiPo) 



,-(logJV(P))V(4t) 



E 

{R}r 



1 



2mR sin Or J^l + e 



re-*'- dr+-tr(/-6(-)) 



-= log 2 - — 



where A„ = i + r,^ are the eigenvalues of A . 

We will now estimate each of the right-hand-side terms of equation l|2.1b for < 
t < 5. 

Parabolic term. 

Lemma 2.1. Let C3 = 8t. Then 



T log 2 r 



— e-'^'Hl + zr) dr + - tr(/ - ©(-)) = 03(7) (0 < i < 5). 



/47rt 27r 



Proof. We have IGR65I page 943 8.361] 

r' 



1 1 



u 1 — e " 



(Re(z) > 0). 



Thus it follows that for ?■ G M, 

1-0(1 + ir)| < |log(l + ir)| + 



du 



|log(l+ir)|- 



1 1 



1 — e " 

where 7 = —0.577... is Euler's constant. Hence 

+ ir) - log(l + ir)\ < -7. 
Now an elementary estimate shows that 
(2.2) |V'(1 + ir)| < 4|r|^/^ +4 (r G M), 

and it follows that 



u 1 - e-« 

du= |log(l+ir)|-(7/.(l)-logl) = |log(l+ir)h7, 



Noting that 



e-*'-'(4|r|i/4 + 4)dr < y (0 < t < 5). 



log 2 1/2 , 
^ < — (0 < i < 5), 



VlTTt t 

and that (5(1/2) is a unitary (t x t) matrix with real entries (hence the diagonal 
entries are ±1), concludes the proof. □ 
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Elliptic term. Next we estimate the elliptic term 



^ 9m „ si 



{R} 



-e-*'' dr, 



where m^j are integers and < Or < t:. 
Lemma 2.2. 



E 



1 



{R} 

In particular, if 



2mR sin 0RJiiil + e-^^^ " f-^ 2mR sin Or \ 29 r 2{ii - Or) 



^ ^ 2mR sin Or ( 20^ + 2(7r - Or) 
iRi 



Then 



E 

{R} 



1 



2mi^sin6li^yR l + e-2'^'■ 



:e-*'■'dr=C'4(-) (0<i<5) 



Proof. For < < tt, and i > 0, 



-2r6l 



1 + e- 



g-2re ^ > g 
g2r(7r-e) ^ ^ 0. 



The lemma now follows from the equation 



2eR 2(7r - Or) ■ 



□ 



Identity term. 

Lemma 2.3. For allt > 

(2.3) 



\T\ 

/ re~^ * tanh(7rr) dr 

47r 



< 



mi 

Att t' 



Proof. Since | tanli(7rr)| < 1 for r e 

(2.4) f re"'''* tanh(7rr) dr = 

47r ./b 



47r 



re * tanh(7rr) dr <2 



47r 



re ^ *■ dr ~ 



47r r 



□ 



Hyperbolic term. Since F is discrete, there exists a hyperbolic conjugacy class 
{Ps\ with minimal norm N{Ps) > 1. Choose c > so that 



1 <c<N{P,). 



Lemma 2.4. For <t < 5, 



1 



E 



logAf(Po) 



/4^rtr ^(P)l/2_iV(P)-l/2 



-(logAr(P))V(4t) < 



745B 
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Proof. First note that for a hyperbolic element P, 

NiP) = NiPoY 

for some integer j > 0. Since N{P) > 1, N{Po) < N{P). Noting that 

log 2: , > 

implies 



1 V- logA^(-Po) 
/4^.m iV(P)i/2_Ar(p)-i/2' 



< 



log N{P) 



/Ant f-^ N{PYn - iV(P)-i/2 



-(logAr(P))V(4t) 



-(log7V(P))V(4t) 



-(logAr(P))V(4t) 



Rewriting the last sum as a Stieltjes integral 



^(logAf(P))V(4t) 



1 



-(logan)V(4t) 



Integrating by parts, and the trivial counting lemma (Lemma II. 3t yield 



/Ant 



-(iog.)V(4t) rf^(^) < 



1 



log a; 



VAnt 
B 1 

B 1 
2t V4^7i 



-(log.)V(4t)^^^^ 



i3('ie*erf(Vi) + ic 



< 



< B e* 



< 



2te 

(logx)e-(i°s^)'/(''*) dx 



/Ant 



/Ant 



745B 



(0 < t < 5). 



The last estimate is a trivial estimate. 



□ 



Spectral terms. We have now bounded all the terms on the right side of the trace 
formula. By Lemmas l2. 1112. 2112. 31 and l2.41 we have 



Lemma 2.5. Let 



Cs = 8t + C4 + — + 7A5B. 
An 



Then 



An 



05{-) (0<t<5). 



3. Explicit bound for the spectral counting function 



We would like to apply Lemma [TTT] to the equation 



(3.1) 



n=0 



An 



_f 2 -(j) 1 
3 — — ( — h ir) ar 
m 2 



0^{-) (0<t<5) 
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but a minor problem stops us: The real valued function —^{\ + ir) may be neg- 
ative for some values of r. However, at the cost of a small error term, we will 
(following Selberg) replace the signed function —^{\ + ir) with a non-negative 
function. 



Notation. Recall from ijl.Sl that Ci, • • • , Cr are a complete list of representatives of 
the r equivalence classes of cusps, and Ai, . . . ,Ar are elements of PSL(2, M) satis- 
fying AjQ = oo. Set 

r = A^rA-\ 

It follows that OD is a cusp for . 

For each 7 G PSL(2, M) let 0(7) denote the lower left-hand matrix entry of 7. Now 
set 

\ = {jer\ 7(c^) ^ 

and let 



(3.2) 



ci = min {|c(7)| | 7 S Ai}. 

1—1...T 



Since F is discrete, Shimizu's Lemma IIEGM981 page 48] asserti3 that ci > 0. 

Modified scattering matrix. Recall from i ] 1.41 that (/)(s) is the determinant of the 
scattering matrix. Let ai > a2 > ■ ■ ■ > o-aa be all the poles of 0(s) in the interval 
(i, 1], and let 



, s - 1 + cr , 



Lemma 3.1. for all r e M, 



(1) -^{\ + ir) is non-negative. 
(2) 



6' 1 . , 
-(2+^0 



< 2|logci| 



Af 



Proof. Item (1) is proved in IIVen82[ Theorem 3.5.5]. (2) follows from elementary 
estimates: Taking the logarithmic derivative of both sides of the equation below 



7 — 1 



plugging in , 



r, and recalling that ^ < cr^ < ■ ■ ■ < ai < I; we obtain 



-h+^r) 



-21ogci+^— T 



|21ogci|+^ 



^ (-2 +^^- ^j)(2 +^^- ^0) 

M 



< 



= 2|logci| 



{^M-W + r^ 



□ 



'ForPSL(2,Z), ci = 1. 
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Let 



W{\) 



0(A) 



0, 
0, 



dr, A > i, 
0< A< i, 



0<A<i 



andletA^(A) ^ |{A„ | A„ < A}|. 
Lemma 3.2. Let 



C7 = 



47r 



2|logci| + 
2|logci| 



V5 

/l67f' 



r/je« the following hold: 
(1) 



d(iV(A) + g(A)) < ^ (0<t<5). 



(2) iV(A) + 0(A) < SCeA (A > i). 

(3) Af(A) + Ty(A) < (3C6 + C7)A (A > i). 



Proof. (1) A simple calculation using A^ 



f-^, shows that 



e-*^d(iV(A) + Q(A)) = e-*/4 
Since e^*/'' < 1, 



\ra=0 



47r 



/.oo / oo 

/ e-*M^(A) + 0(A))<^< 



An 



-"^-^(|+-) 



— ( — \- tr) ar 

r 2 , 



Note that both sides of the above inequality are positive. Let 
(3.3) 

Then rewriting | |3.1| I gives us 



(3.4) 



47r 



-^{^-+^r) + U{r)] dr 



<9l (0 < i < 5). 



But Lemma [3.11 and a crude estimate (see equation | |3.6|| for the idea of the esti- 
mate), show that 



47r 



— / e-*'- \U{T)\dr< 2|logci| + 



1 1 



(o-AT - 5)V Vi&^Vt 



and since 
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it follows that 
(3.5) 



< 



C5 + (2|logci| 



t 



(0 < t < 5). 



(1) is now proved. 

(2) follows from that fact that N{\) + Q{\) is non-decreasing, (1), and Lemma [TT 
To prove (3), note that | |3.3t implies that 



N{\) + W{\) < N{X) + Q{X) 



1 

47r 



Next, using Lemma [3. 11 a trivial estimate gives 



(3.6) 



\U{t)\ dt < 



2|logci| 



\U{t)\ dt. 



< 



2|logci| 



dt 
N 



dt. 



(3) now follows from the equation 



A + -<4A (A>l/4). 



3.1. Implied constant for the spectral counting function. Recall that A„ 
Also note that for A„ > |, we have r„ > 0. 

Let 

1 



□ 

„2 



(3.7) 



/i(r) ^ |{r„ I < r„ < r}\ + 



An 



dt. 



Lemma 13^ 3) immediately implies that 



^l{r) < N{r' + 1) + Wir^ + i) < (SCg + Cj)ir'- + 1) (r > 0). 
We summarize our result: 
Theorem 3.3. Let C = SCe + C7. Then 



^i{r) <C^' + -j (r > 0). 

3.2. Explicit estimate for the modular group. Let F = PSL(2,Z). In this sec- 
tion we will give an explicit upper bound for the number C which depends on 

C3, C4, C5, Cg, and C7. 



The constant B. From the trivial counting lemma (fjLS) we have 

min(2.r) 

B = A'kY{Y + 1) ^ e^''^ 

where r = 1 . Let F be the standard fundamental domain of F 

J" = {z G H I - ^ < Rcz < i, |z| > 1}. 
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The euclidian distance of to the x-axis, is e = We can decompose T 
J"oUJ'(r), where 



Tiy) = {z G H I - ^ < Rez < ^, Imz > F}, 



with Y ^2. 



Recall that A'l = {z e J" | | < Im(z) < y + 1}. Let be the hyperbolic diameter 

of K\. To estimate di note that cosh(p(z, w)) = 1 + jj^t^^j-j;, where w) is the 
hyperbolic distance between z and w. An elementary calculation shows that d\ < 
2.3. Hence B < 753. 

The constant C4. It is known that F has two classes of elliptic elements, repre- 
sented by 

'0-1 
1 



R = 
and 



S 



1 -1 
1 

Elementary calculations show that nij^ = 2 (the order of the centralizer of R in F) 
and that = 2cos(6'7^) where < Or < tt; thus 9ji = ^. Similarly ms = 3 and 
6*5 = |. From Lemma IZ2l it follows that C4 < |. 

The remaining constants. Now C3 = 8t = 8, and | J^l = |. 

The number ci = 1 and there is only one pole (at s = 1) for the scattering matrix 
0(s), so A/" = 1 and a/^ = 1. Now CqjCj, and C can be easily computed. The result 
is 

Theorem 3.4. Let F = PSL(2, Z). Let ^{r) be the spectral counting function I I3.7I I, and 
let 

C= 1,682,997. 

Then for r > 0, 

■ 2 1 



The estimate for C is probably not sharp. However, we did not use any of the 
number theoretic properties of PSL(2,Z) to arrive at the estimate; we only as- 
sumed that F was a cofinite Fuchsian group. 

4. Bound for the Huber constant 

In this section, we adapt Randol's proof ( IRan84l ) of the prime geodesic theorem 
to an arbitrary cofinite Fuchsian group. Along the way, using Theorem 13.31 we 
will obtain an upper bound for the Huber constant Cm . Our main references are 
||Ran84[ pages 295-300] and |Bus92[ Section 9.6]. 



4.1. Explicit bump function. For T > set 

It{x) -- 



1 \x\<T, 
\x\ > T. 



For a functions /, g on M, we set 



fix) ^ -= I f{y)e-^^y dy 
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and 

{f*9){x) = -^ I f{x~y)g{y)dy. 



/2tt 

Lemma 4.1. There exists a nonnegative function (j){x) with the following properties: 

(1) e C(f (K) with supp(0) c [-1, 1]. 

(2) m = 1- 

(3) Let Ci = 2(e - 1). Then 0(ez) - 1 = Oi(e)for 0<e<l,andze 1]. 

(4) Let b>0, = Then for z = r + it, \t\ < b, we have 

0(z) = O2((l + |rr2), (r e K, |i| < 5). 

Proof. Let 



Co = y cxp \^— — -j dx = 0.4439938 . . . , 

and let 

^ / exp(^) \x\ < 1, 



^^"^ CO 1 \x\>l. 
(UJ and (|2]l now follow from the definition of (p. 

For w e C, |w| < 1, observe that 

w'^ w'^ \w\ IwP 1 1 

|e--l| = \nj+ — +—+. . . , I = |^„|.|1+U+ M_^. . . , I < + . . , I = 

For < e < 1, and z e «[— 5, 1], note that |ez| < 1. Let w ~ ez. By definition, and 
©, it follows that 

- 1 = ^ / 0(y)e-™« - 1 = 
V 27r JR 

^ ^ ^{V)e-^^y dy--^ I ct>{y) dy = C <^(y)(e-™^ - 1) dy. 



\/27r Jh v27r Jr v27r 

Now, since 



1 



#<i (.eM), 



it follows that 



\4>{tz)-l\< j |e-™^-l|dy< j (e-l)|u;|d2/ = 2(e-l)|u;|=2(e-l)|ez|<2(e-l)e. 
This proves 

Next we prove iSll. Suppose | Im(z)| < h. Then since (/)(0) = 1 it follows that 

\m\<e\ 

Integrating by parts, twice, shows that for 2; 7^ 0, 
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Now, a simple calculatiorQ shows that 

|0"(2/)|<1O6 i~l<y<l). 

Hence, for z ^ 0, 
Finally, observing that 



and 



it follows that 



1 2 



,M 1 424 . 1 424 . 1 848 >, 

- + 1 V2^ -r2 + lV2^ -(|r| + 1)2^2^ 
since ^ ^ 

T^TT - (|r| + l)2 

This concludes the proof. □ 

4.2. Explicit lemmas. In this section we prove some lemmas which allow us to ob- 
tain explicit constants from Randol's proof of the prime geodesic theorem. Recall 
(SiU/ and the definitions of It{x) and cj). For e > 0, define 



and define 
Set 

It follows that 
where 

Note that 



g'j.{x) = (2cosh(|)) (/T*0e)(x). 



S{w) = I - sin (Tw) ) (j)e{w), S{0) = 2T. 



(t)^{w) = cj){ew) 



The main idea behind the proof is to apply the Selberg trace formula to the pair 
h^{r), g^ix)', and to extract the leading and error term from the trace formula. 
Much of our effort will be spent on giving explicit estimates to each term in the 
trace formula. 

We next set up some important notation. For < Afc < -j , a small eigenvalue, define 
rfc, and Sfe, by 

Tk = « - Afe = i{sk - i). 



%irst note that 4>"i^) = j-^t!f j)4 (3^'* — !)■ Next using elementary calculus maximize ^^^4 
and (3x'* — 1), separately, and multiply their maximums. Of course, this argument does not yield the 
optimal bound. 
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Let A denote the number small eigenvalues: 

A^\{Xk\Xk<\}\. 
Explicitly, = Aq < Ai < • • • < < i. 
Lemma 4.2. Let do = 10-^ V^, then for r > 0, < e < 1, 



(1) \hf,{r)\<C,oe^/Hl+r)-\l + er)-\ 

(2) Forre [0,i], 

T/2 

\hfrir)\ < Ciof— 
1 + r 

(3) For r e [ioo), 



\hUr)\<Cio 



eT/2 



Proof. We prove (1); (2) and (3) will then follow from elementary estimates. 
Recall that for r > 0, 

hUr)=Sir+'-) + Sir-'-), 

where 

Slw) = [ — sin(Tw) I (j)f(w), 
\w J 

and note that 0e(w) = (j){ew). By Part|4]of Lemma l4Al with b — 1/2, it follows that 

?(6(r±^))==02((l + er)-2), 



where 



Now 



and 



Hence, 



V 27r 



2 5 

< 



|r±§| - 1 + r' 



|sinr(r±^)| <e^/^ 



|/i|^(r)|<2-^C2(l + 6rr2. 
1 + r 



□ 



The next lemma will be used repeatedly. 

Lemma 4.3. Let u{r) be a non-decreasing function on [0, oo) that satisfies 
\u{r)\ <cr^ + b (r > 0); c,6 > 0. 

Then 

/ \h'j,{r)\ duir) < K 

Jo £ 

where K = do (6c + i§^) . 
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Proof. By Lemma [4.2[ it suffices to estimate 



(4.1) 



/ du{r) and / — du{r). 

Jo 1 + r Jiu T 



Integrating by parts, yields 



' Mr) 



1 + r 



1 + r 



u(r) 



c/e2 + b 



+ b + 



1 + 1/e 

A similar argument shows that 

poo -j^ 



{i + ry 

^1' cr^ + b 







(l + r)2 



dr < 



c , 3c b 4c + 26 
dr < - + b-\ h-< . 



3 duir) < (c + b)e + (c + = (2c + ^)e. 

l/e r 6 6 



Hence 



/■°° / 106 

\hfrir)\duir)<C,o[6c+ — 



106 \ 



□ 



4.3. The leading term of the prime geodesic theorem. In this section we show 
that the sum 

gives us the leading term of the prime geodesic theorem. In the next section, we 
will see that the above sum is a term in the Selberg trace formula. 

Lemma 4.4. Let e = e^-^/'*, and 

Ci2 = (^-1) fl + 3Ci 



1-si 



Then 



where 



Sk 



(1 + Ci) + 2Ci + 2. 



Oi2(e3^/^), 



Sl 



i - Al Al < i, 



Ai > i 



Proof. First we deal with the case Aq = 0. We have 

hriro) = h'ri^) = - sin(ir)0(ie)+2T = ^(ie)(e^-e-^)+2T = (l+Oi(e))(e^-e-^)+2r. 
Hence, 

|/i^(ro)-e^| < Ciee^ + e-'^\Oi{e)\+2T <Ciee^ + l + Cie + 2T < (2Ci + 2)e3^/^ 

The last inequality follows from the trivial bound 1 + 2T < 2e'^^/^, and by noting 
thate = e-'^/^. 

Next we handle the (possible) case of nonzero small eigenvalues. Let < fc < 

A = Afc, r^rk, s = Sk- 
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It follows that 



' s — 1 



Here, ^ < s < si <l.By Lemma OJ Part|3l 



Therefore, 



hUr) 

s 



(ies) = 1 + 01(6), and 0(ie(s - 1)) = 1 + Oi(e). 



< + (1 + I 0,ie)\)e-^^ + ^—(1 + + | Oi(e)|) < 

s \s-l\ 

\e^Cie + (1 + Ck) + —^{l + e^/2)(l + Cie). 
k 1 - si 



Next, set e = e We obtain 



Then, finally, we have 



< 



l + 3Ci + -^(l + Ci) e3^/4 
1 - Sl / 



< 



{A -I) l + 3Ci 



1-si 



(1 + Ci) +2C1 + 2 



Of course, if there are no nonzero small eigenvalues, A = l. 



p3T/4 



□ 



4.4. Application of the trace formula. Let F be a cofinite Fuchsian group, and let 
M = F \ H be the corresponding hyperbolic-orbifold. Let C{M) denote the set of 
closed geodesies of M, and let V{M) denote the set of prime (or primitive) closed 
geodesies (see IBus92[ page 245]). For each 7 e C{M) there exists a unique prime 
geodesic 70 and a imique exponent m > 1 so that 7 = 7™ . Let / (7) denote the length 
of 7. Associated to 7 is a imique hyperbolic conjugacy class {P-y^r with norm 



For 7 e C{M) define 



A(7) = log7V^„ =Z(7o). 



Let 



7eC(M) ^^7 ^^7 



H{T)= E A(^)y 



1 + iV-i 



;(7)<T 



Note that H^{T) is an approximation of H{T), and that for any e > 0, ( l|Ran841 
page 298]) 



(4.2) 



H,{T~e) < H{T) < H,{T + e). 
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An application of the Selberg trace formula yields 
(4.3) 



47r 7r (j) 2 47r 7^ 



+ i/,(r)+^- f \ hf,{r)dr+'-^triI -&{-)) 

f-f 2miiS\neR 1 + e ^'^'^ 



{i?,} 



/i|,(0) 



rff^(O) log 2 - ^ / /iT + dr. 



Our next goal is to let e = e ^, and estimate each term, in terms of the variable T. 
Recall from Theorem |3.3| that 



^i{r) < C + - (r > 0) 



Lemma 4.5. Let e = e 



Ci3 — -^CCio- 
6 



Then 



Proof. 



47r7R 2 



< 



dr = / |/i|,(r)| dix{r). 



By Theorem|33l ^^(r) < C{r'^ + \). The result now follows from Lemma [4.31 □ 
Lemma 4.6. Let e = e"^/"', 



2Q6r T 
C^i4 = Cio^ + Cio- + 2Tlog2, 



then 



^ I \hfr{rmi+ir)\dr + 



^..■(/-s(i)) 



Proof. We previously saw (equation \2.2\i ) that for r e M, 

\^{l+ir)\ < 4|r|^/^ + 4. 
Hence \ip{l +ir)\< 8(1 + \r\), and 

/.OO /-OO /-OO 

— y_ |/i^(r)V'(l + ^r)\ dr < — j _ |/i^(r)|8(l + |r|) dr ^ — |/i^(r)| rf«(r), 
where 

/"■ 3 
u{r)= 8(1 + \x\) dx = 16{— + r) < 16{-r^ + 1) (r e M). 

J — r ^ ^ 
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An application of Lemma l4!3l yields 



T 

2^ 



\h'T{r)\ du{r) < do 



296r e^/^ 



Next, an application of Lemma l4!2l the inequalities |tr(/ — 6(i))| < 2r and 
135,(0)1 < 2, imply 



^M/-e(i)) 



|r5!,(0)log2| < I§^Ie^/2+2rlog2 < 



2rlog2)e3^/^ 
□ 



Lemma 4.7. Lef e = e''^/^, 



^15 



56Cio 



E 

{R} 



1 



2to /? sin 9fi 



then 



E 

{i?} 



1 



2mB^ sin 9r J^l + e 



-2r0f 



— hfp(r) dr 



= Oi5(e3^/^). 



Proof. Since 



E 

{R} 



1 



2mfl sin Jr 1 + e 



-2ref 



l + e- 



:h^{r) dr 



< 



< 1, 



E 

{R.} 



1 



2niR sin 



|/i|,(r)| dM(r), 



where u{r) ~ 2r. Hence \u{r)\ < 2r^ + 2, and applying Lemma l43l yields 



E 

{R} 



1 



2m_R sin 6'_r Jr 1 + e 



:h^{r) dr 





f 56Cio 







E 

{R} 



1 



2mfl sin 6*7? 



p3T/4 



□ 



Lemma 4.8. 



^ / |/i|,(r)rtanh(^r)| dr < -^|J-|Cio . 

Att 27r e 



Proof. Since r tanh(7rr) is even, it suffices to estimate 

LIJ^ \hUr)rte.nHnr)\dr<Ll \hf^^r)\r dr = LI \hUr)\du{r), 

where u{r) = x dx = ^ . (Note that we used the trivial estimate | tanh(7rr)| < 1.) 



The result now follows from Lemma [431 

Our main result for this section is 
Lemma 4.9. Let e = e~-^/^, and let 



□ 



Cie ~ C12 + Ci3 + Ci4 + Ci5 + 7^— |-7-'|Cio- 

ZTT 
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Then 



Proof. The proof follows immediately from equation (|4.3l l, and Lemmas 14.41 14.51 

□ 



4.5. Upper bound for the Huber constant. In this section we prove the main re- 
sult of this paper. 

4.6. Asymptotics for H{T). 
Lemma 4.10. Let e = er^l^, 

Ci8 =4^ + 5Ci6. 



Then 



(1) H{T)-H,{T)^On{e^T/^). 
(2) 



Proof. We first prove (1). Recall from equation l|4.2b that 

(4.4) H,{T ~ e) < H(T) < H,(T + e). 
In addition, it foUows from the definition of (T) that 

(4.5) H,{T - e) < H,{T) < H,{T + e). 
By Lemma [Ol we have 

H,{T±e)^ V ^ + Oi6(e3™/4)^ 

— Sk 

Now, setting e = e~^/^, it follows from elementary estimates that 

\H,iT + e) - H,iT - e)| < iAe^^/* + ^Cx^e-'^l\ 

Part (1) now follows from (|4.4|l and l l4.5l l. Part (2) follows immediately from Lemma l49l 

□ 

4.7. The Chebyshev function. Define 

^{T)^ A(7), 
K7)<r 

where A(7) = Z(7o). Since F is discrete, there exists a hyperbolic conjugacy class 
{Pg} with minimal norm N{Ps) > 1. Choose c > so that 

1< c < N{P,). 
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Lemma 4.11. Let 



Cig — Ci8 



^i2A + Cis). 

I — c ^ 



Then 



— Sh 



Proof. By Lemma |4T0l 



^ Sfc 



Since 



we can rewrite H{T) as 



lij)<T ^ 



i?(r) = vi;(r) + ^ A(7)-^. 

Both of the above terms are positive for all T. Thus, from Lemma l4.10l and a crude 
estimate (for each small eigenvalue, ^ < < 1, hence e'''^* < e*.), it follows that 



(4.6) 
But 



2N~'^ 9 9 



;(7)<T 



—{2A + Cis){l+T). 

1 — c ^ 



The last inequality is obtained by integrating by parts and using | |4.6t . Next, using 
the crude boimd (1 + T) < 2e^'^/*, we obtain 



2N-^ 



ih)<T 



Finally 



*(r)-E- 



aTSk 



Sk 



< 



(Ci8 + - -i2A + C^s))e'^/'- 

1 — c ^ 



□ 



4.8. The Huber constant. Let 

e(r)= 

ih)<T 

7 primitive 

Choose fi> Oso that < ^(7) for all 7 G C(A/). 



AN EFFECTIVE BOUND FOR THE HUBER CONSTANT FOR COFINITE FUCHSIAN GROUPS 



Lemma 4.12. 



where 



^ — ^ J. 



C20 — CiQ H — {2 A + Cis). 
A' 

Proof. Let m(T) = [T/n], then ( IIBus92[ page 251]) 

m(T) 



T 



*(r) - e(r) = V err/m) < m(T)e{T/2) < -*(t/2). 



By sal. 



-*(T/2) < -Te^/^{2A + Cis) < -e^^^H'^-A + Cis). 

Hence 

|e(r) - *(T)| < -e^^/^{2A + Cis). 
The lemma now follows. 



Let X ~ , and let 



9{x)^ A(7). 



7 primitive 



With this change of variable, Lemma l4.12l becomes 



Let 



li(x) 



dy 

logy' 



and recall the definition of the constant c (i l4.71 . We have 
Lemma 4.13. Let 



1^ 
'log 2 



Then for s e [|, 1] 



/c logy 

Proof An elementary argument shows that 



logy 



dy - \i{xn 



|2- 


logc 




< C21. 




dy 




Jc^ logy 



Now, if > 2, then 



If < 2, then 



dy 



logy 



- ^ ^ log 2 - ^ ^ log 2 



dy 



logy 



< (2-c«)^ < (2-ci/2)J_. 
~ ^ ^logc" ~ ^ ^logc 
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Hence, in either case. 



dy 



logy 



<\c-2\-l- + \2-c'/^\-^. 
log 2 log c 



Lemma 4.14. Let 



Then for x > 1, 



C.o= 1- 



log2 



logx 



□ 



Proof. Integrating by parts give us 

ii(^)^ r -^<-j^+ r ^ 



1 



^ ^ dy 

1 2 logy ~ loga- ' J 2 (logi/)2 ~ logi- ^ log 2^2 logy' 

The lemma now follows. 
Theorem 4.15. Let 

c3/4 3 

Cu = C21A + C20- \- C20 + C2aC22 + 7(^20(^21 • 

log c 4 



Then for all x > 1, 

In particular, Cm < C*!/ • 
Proof. First note that 



logx 



7r(x) 



logy ■ 



First, for notational convenience, set 



f{y) EE e{y) ^Y.y—= 02o(y3/4). 

Afc<T 



Next, 



By Lemma |4l3l 



Now 



logy logy 



E 



logy 



, ,^ logy 

" dfjy) ^ fix) f{c) 
log y log X log c 



/(y) 



dy 



y(logy)2 



□ 
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and (using Lemma [4.13|l , 



\fiy)\ 



dy 



< C20 / y 



„,3/4 dy 



y{\ogyY 



C2~ 



,3/4 ^3/4 3 



log c log X 4 



y 



1) dy 



logy 



But by Lemma I4l4l 



■ li(a;3/4) < C: 



„3/4 3/4 o 

lege log a; 4 

^3/4 



4 ^ ' ^ logo; 
Hence, putting everything together, for a: > 1 , 



^3/4 ^3/4 
< C21A + C2Q- YC20-. 



log c log X 



a;3/4 3 

■C20C22-. h 7(^20^21 • 

log a; 4 



And using the trivial inequahty 1 < ^j—^ , we obtain 



logx ' 



where 



Cu — C21A + C; 



20 



c3/4 3 

1 \- C20 + ^^20^22 + -C2oC21- 

log c 4 



□ 



4.9. The Ruber constant for PSL(2, Z). For T = PSL(2, Z), we can take c = 6.85. 
Also, r has no small eigenvalues (except for the trivial one), so ^ = 1 . Using the 
computations in i )3.2[ an explicit value for Cu can be calculated: 

Theorem 4.16. Let V = PSL(2, Z). Then Cu = 16,607,349,020,658 is an upper hound 
for the Huber constant. 

Using the theory laid out in |Sar82L and the excellent open-source computational 
program PARl/GP IITheOSI , we can explicitly calculate (and list out) the length 
spectrum for PSL(2, Z). Our computations suggest that 

Cm < 2. 

So it seems that our result is probably not sharp! 



Computer experiments for PSL(2, Z) 



Let Fo ~ PSL(2, Z), and F be a finite index subgroup of Fq. Let {Po} denote the 
set of primitive hyperbolic conjugacy classes of F. We compute length spectrum 
(norms of primitive h5^erbolic conjugacy classes) using the ideas of IISar82ll . 

Let Q{x,y) = ax"^ + hxy + cy"^ be a primitive, indefinite binary quadratic form of 
discriminant d = iP' — Aac > 0. Two forms [a, 6, c] and [a', h' , c'] are called equivalent 
(in the narrow sense) if they are related by a imimodular transformation. 

Consider the map which sends the form [a, 6, c] to the matrix 
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where td, > is the fundamental solutions of Pell's equation t"^ — du^ = 4. Set 

ed = ^ . 

Let V be the set of fundamental discriminants {d > 0,c? = 0, 1 mod 4} and let 
h{d) denote the number of inequivalent primitive forms of discriminant d. For each 
d e I? let Qi (d) , . . . Qh{d) (d) be a complete set of inequivalent forms of discriminate 
d. Finally, set 

QF= \J{Qiid),...Qhid){d)} 

dev 

Lemma 4.17. f lSar82l ) The set (t){QF) is the set ofconjugacy classes of primitive hyper- 
bolic elements o/PSL(2, Z). The norm of each hyperbolic element (j){Qi{d)) is e^. 

Using the above ideas, one can compute the primitive conjugacy classes for Fq = 
PSL(2, Z). We next explain how to compute the length spectrum for F, a finite- 
index normal subgroup of Fq . 

Let F < Fq be a finite-index normal subgroup of Fq = PSL(2, Z), G = Fq/F, and 
n = \G\. For each primitive conjugacy class {P^} of Fo let m represent the order of 
Py e G. It follows that {PI^} is a subset of F that splits into conjugacy classes with 
respect to F. 

Lemma 4.18. f lSar82| ) The number ofconjugacy classes into which {P™} splits in F is 

k = n/m = \G\/m. 

5. Algorithm for computing primitive conjugacy classes 

Let F O Fq, X > 0. We wish to compute the norms and multiplicities of all primitive 
h5rperbolic conjugacy classes {P-y} with norm N{{Py}) < x. 

Step 1 . Compute the primitive conjugacy classes for the modular group Fq of norm 
X or less. Since is the norm of a prime element that is mapped by the function cj) 
it follows that ^/x > Cd > Vd/2. Thus it suffices to consider discriminants d < 4a; 
in order to guarantee that all norms less than x are achieved. In other words, if 
d > Ax then > x. 

For each d £ = {d £ V \ e^^ < x} , compute h{d) and representatives of inequiv- 
alent quadratic forms of discriminant d. For each representative Qi{d) compute 
(f>{Qi{d)) in order to obtain representatives for {P-^} in Fq. 

Step 2. From step 1, let z = P~f he a representative of a class in Fo of norm e^. 
Compute the order m of z in G = Fq/F. Then {P™} splits into k = |G| /to primitive 
conjugacy classes each with norm e^"' . 

Step 3. From step 2, we count the primitive conjugacy classes with norm e^™ < x. 

One thing we can do is verify the prime geodesic theorem for the group F. We can 
implement these algorithms using the computer package FARI/GP. 

Pari/GP Algorithm. 

: IsProp{D) ^if{{D%A^=l\\{D)%A==0) && issquare{D), 0,1); 
: QFBclassno{D) = qfbclassno{D) *if{D < || nor'm{quadunit{D)) < 
0,1,2); 
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QUADunit{D) = quadunit{D) * if{D < || norm{quadunit{D)) < 
0, quadunit{D), 1); 
Li{x) = —eintl{log{l/x)); 
MaxD = 10000; 
{ 

for{D = 5, MaxD, if{IsProp{D) == 0, 

n = QUADunit{D) A 2.0; 

if{n < MaxD, 

m = QFBclassno{D); 

q = Li{n); 

writeC FILE. CSV" , n" , " m" , " D" , " q)))); 
} 
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Table 4. Data for the modular group 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


6.85 


1 


5 


4.68 


1 


4.68 


3.68 


3.05 


1.21 


2.62 


1.41 


13.93 


2 


12 


7.75 


3 


2.58 


4.75 


4.44 


1.07 


3.73 


1.27 


22.96 


2 


21 


10.87 


5 


2.17 


5.87 


5.92 


.99 


4.79 


1.23 


33.97 


2 


32 


14.17 


7 


2.02 


7.17 


7.49 


.96 


5.83 


1.23 


33.97 


1 


8 


14.17 


8 


1.77 


6.17 


7.49 


.82 


5.83 


1.06 


46.98 


2 


45 


17.69 


10 


1.77 


7.69 


9.15 


.84 


6.85 


1.12 


61.98 


4 


60 


21.45 


14 


1.53 


7.45 


10.87 


.68 


7.87 


0.95 


78.99 


2 


77 


25.45 


16 


1.59 


9.45 


12.68 


.75 


8.89 


1.06 


97.99 


2 


24 


29.69 


18 


1.65 


11.69 


14.55 


.80 


9.9 


1.18 


97.99 


4 


96 


29.69 


22 


1.35 


7.69 


14.55 


.53 


9.9 


0.78 


118.99 


1 


13 


34.17 


23 


1.49 


11.17 


16.48 


.68 


10.91 


1.02 


118.99 


2 


117 


34.17 


25 


1.37 


9.17 


16.48 


.56 


10.91 


0.84 


141.99 


4 


140 


38.9 


29 


1.34 


9.9 


18.48 


.54 


11.92 


0.83 


166.99 


4 


165 


43.86 


33 


1.33 


10.86 


20.53 


.53 


12.92 


0.84 


193.99 


2 


48 


49.06 


35 


1.40 


14.06 


22.65 


.62 


13.93 


1.01 


193.99 


4 


192 


49.06 


39 


1.26 


10.06 


22.65 


.44 


13.93 


0.72 


223 


4 


221 


54.49 


43 


1.27 


11.49 


24.82 


.46 


14.93 


0.77 


254 


2 


28 


60.15 


45 


1.34 


15.15 


27.04 


.56 


15.94 


0.95 


254 


4 


252 


60.15 


49 


1.23 


11.15 


27.04 


.41 


15.94 


0.7 


287 


4 


285 


66.04 


53 


1.25 


13.04 


29.31 


.45 


16.94 


0.77 


322 


2 


80 


72.17 


55 


1.31 


17.17 


31.63 


.54 


17.94 


0.96 


322 


1 


20 


72.17 


56 


1.29 


16.17 


31.63 


.51 


17.94 


0.9 


322 


4 


320 


72.17 


60 


1.20 


12.17 


31.63 


.38 


17.94 


0.68 


359 


4 


357 


78.51 


64 


1.23 


14.51 


34 


.43 


18.95 


0.77 


398 


2 


44 


85.08 


66 


1.29 


19.08 


36.42 


.52 


19.95 


0.96 


398 


8 


396 


85.08 


74 


1.15 


11.08 


36.42 


.30 


19.95 


0.56 


439 


2 


437 


91.88 


76 


1.21 


15.88 


38.88 


.41 


20.95 


0.76 


482 


4 


120 


98.89 


80 


1.24 


18.89 


41.39 


.46 


21.95 


0.86 


482 


8 


480 


98.89 


88 


1.12 


10.89 


41.39 


.26 


21.95 


0.5 


527 


4 


525 


106.12 


92 


1.15 


14.12 


43.94 


.32 


22.96 


0.62 


574 


4 


572 


113.57 


96 


1.18 


17.57 


46.53 


.38 


23.96 


0.73 


623 


2 


69 


121.23 


98 


1.24 


23.23 


49.16 


.47 


24.96 


0.93 



f. Length spectrum calculations for the modular group. 

(1) X, the actual primitive length spectrum (norms); when norms repeat, they 
from different quadratic forms 

(2) multiplicity of given norm 

(3) discriminant of corresponding quadratic form 

(4) li(a-), the integral logarithm of x where x is equal to the value of column 1 

(5) ■k{x), a running total of the number of norms less than x where x is equal to 
value of column 1 

(6) ; this column verifies the prime geodesic theorem 

(7) \\\lx)-^{x)\ 

(8) x^/^/{\o^(x)Y/^ 

(9) the ratio of (7)/(8) 

(10) x^/^, this is the conjectured true error term 

(11) the ratio of (7)/(10) 
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